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Abstract 

We consider the adiabatic regime of two parameters evolution semigroups generated by 
linear operators that are analytic in time and satisfy the following gap condition for all times: 
the spectrum of the generator consists in finitely many isolated eigenvalues of finite algebraic 
multiplicity, away from the rest of the spectrum. The restriction of the generator to the spectral 
subspace corresponding to the distinguished eigenvalues is not assumed to be diagonalizable. 

The presence of eigenilpotents in the spectral decomposition of the generator forbids the 
evolution to follow the instantaneous eigenprojectors of the generator in the adiabatic limit. 
Making use of superadiabatic renormalization, we construct a different set of time-dependent 
projectors, close to the instantaneous eigeprojectors of the generator in the adiabatic limit, 
and an approximation of the evolution semigroup which intertwines exactly between the values 
of these projectors at the initial and final times. Hence, the evolution semigroup follows the 
constructed set of projectors in the adiabatic regime, modulo error terms we control. 

1 Introduction 

Singular perturbations of differential equations play an important role in various areas of math- 
ematics and mathematical physics. Such perturbations typically appear when one considers 
problems that display several different time and/or length scales. In particular, the semiclas- 
sical analysis of quantum phenomena and the study of evolution equations in the adiabatic 
regime lead to singularly perturbed linear differential equations which are the object of many 
recent works. See for example the monographs Jl], ^31) [201; IIHI- The description of cer- 
tain non conservative phenomena with distinct time scales also gives rise to non-autonoumous 
linear evolution equations, which are more general than those stemming from conservative 
systems, and whose adiabatic regime is of physical relevance, see e.g. [^, [^, jU, [H^ . 

m, la, 0, m. 

The present paper is devoted to the study of general linear evolution equations in the 
adiabatic limit under some mild spectral conditions on the generator. The chosen set up is 
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sufficiently general to cover most applications where the time dependent generator is charac- 
terized by a gap condition on its spectrum. Let us describe informally our result, the precise 
Theorem being formulated in Section |^ below. 

We consider a general linear evolution equation in a Banach space B of the form 

i£dtU{t,s) = H{t)U{t,s), U{s,s)=I, s < t e [0, 1] (1.1) 

in the adiabatic limit e ^ 0^, for a time-dependent generator H(t). This equation describes 
a rescaled non-autonomous evolution generated by a slowly varying linear operator H{t). The 
evolution operator U{t, s) evidently depends on e, even though this is not emphasized in the 
notation. 

The generator H{t) is assumed to depend analytically on time and to have for any fixed t a 
spectrum a{H{t)) divided into two disjoint parts, a{H{t)) = cr(t)Ucro(t), where a{t) consists in 
a finite number of complex eigenvalues a{t) = {Ai(t), A2(i), • • • , A„(t)} which remain isolated 
from one another as t varies in [0, 1]. Moreover, the spectral projector of H{t) associated with 
o"(t), denoted by -P(t), is assumed to be finite dimensional. The part of H{t) which corresponds 
to the spectral projector Po{t) associated with cJo(i) can be unbounded, bounded or zero. In 
the first case we need to assume H{t) generates a bona fide evolution operator. 

This spectral assumption, or gap condition, is familiar in the quantum adiabatic context 
where i3 is a Hilbert space on which H{t) is further assumed to be self-adjoint, see [101, ^5] . 
innji fo'^ example. Note that it is still possible to study the quantum adiabatic limit by 

altering the gap condition in different ways, as shown in [S], |TH], ^3]) IMj; [Bli HI- 

By contrast to previous studies of similar general problems J2], [22j, |2H], [22], we do 
not assume that the restriction of H{t) to the spectral subspace P{t)B is diagonalizable. Such 
situations take place in the study of open quantum systems by means of phenomenological 
time-dependent master equations, [SS], [SHI, W^i EH- We come back to the approach of j35j 
below. 

Therefore, for the part H{t)P{t) of the generator, we have a complete spectral decomposi- 
tion 

n 

H{t)P{t) = Y,Ht)PAt) + Dj{t), (1.2) 
i=i 

where the Pj{t)^s are eigenprojectors and the Dj(tys are eigenilpotents associated to the eigen- 
value \j{t) that satisfy 

n 

Y,Pjit) = Pit)^ Pj{t)Pk{t) = S,kPj{t), and D,{t) = Pj{t)Dj{t)Pj{t). (1.3) 
i=i 

In case ;S is a Hilbert space on which H{t) is self-adjoint or if H{t) is diagonalizable with real 
simple isolated eigenvalues only, the evolution U{t, s) follows the instantaneous eigenprojectors 
Pj{t) in the adiabatic regime in the sense that 

C/(t, s)Pj{s) = Pj{t)U{t, s) + 0(e), as e ^ 0, (1.4) 

as shown in ^01 , 113] i IHOji Qi and ^2j, [2Hj, [221, for example. In other words, transitions 
between different spectral subspaces are suppressed as e ^ 0. This relation remains true 
for certain eigenprojectors if the eigenvalues are allowed to have negative imaginary parts, 
|,S2j . This fact is also well-known and crucial in the study of the Stokes phenomenon 
appearing in singularly perturbed differential equations ^A^■. under analyticity assumptions, 
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one considers certain paths in the complex t-plane, called canonical of dissipative paths, along 
which an equivalent of (|1.4j) is true in order to get bounds on, or to compute exponentially 
small quantities in 1/e stemming from singularities in the complex t-plane. Such methods are 
used in [20], |2j, [SH] and JH], to bound or to compute exponentially small transitions in the 
adiabatic limit when the relevant eigenvalues are real on the real axis. 

However, when eigenilpotent are present in the spectral decomposition 1)1. 2|) . the relation 
()1.4)) cannot hold in general, even for real valued eigenvalues Xj{t). Indeed, the transitions 
between spectral subspaces Pj{t)U{t,0)Pk{0) are typically exponentially increasing as e — > 0, 
rather than vanishing as e. An explicit example of this fact is provided at the end of the 
Introduction. We come back to this mechanism below. 

In this context, our main result reads as follows. We construct a different set of time- 
dependent projectors Pj ^^\t) which approximates the eigenprojectors Pj{t) in the adiabatic 
regime e — > 0. And we show that the evolution U{t,s) can be approximated by a simpler 
evolution, which exactly follows the constructed approximations Pj ^'^\t) of the 

instantaneous eigenprojectors. In other words, we restore the expected adiabatic behaviour 
by trading the instantaneous eigenprojectors for other nearby projectors in the limit e — > 0. 
Note that since the eigenvalues need not be real in general, we also have to take into account 
the contributions stemming from the "dynamical phases" e~*-'^s >^i{u)du/£ ^j^j^]^ gg^j^ ]-,g expo- 
nentially increasing or decreasing as e ^ 0. In case H{t) is unbounded, we assume the part 
H{t)PQ{t) generates a semigroup bounded by \e~^^i> \o(u)du/e]^^ some function Ao(t). 

More precisely, for all j = 1, • • • , n and for any < t < 1, we construct perturbatively a set 
of projectors close to the spectral projectors of H{t), see Sectional 

n 

Pp'\t) = P,{t) + 0{e), and pf'^'\t) = I - (1.5) 

i=i 

Let W'^*^'^\t) be the intertwining operator naturally associated with the projectors ^^\'t)j 
k = 0, • • • ,n introduced by Kato jlH], such that 

W''*^'Ht)Pl*'^'\0) = pf^'\t)W'^*^'\t), k = 0,---,n. (1.6) 

The approximation is then defined by 

(t, 0) = 1^«*(") (t)^"^^'^ {t, 0) (1.7) 

where <I>^*(^^(t, s) commutes with all the P^ ^^\o) , k = 0, ■ ■ ■ , n for any t and satisfies a certain 
singularly perturbed linear differential equation, see (|6.17)) below, which describes the effective 
evolution within the fixed subspaces P^ ^^\o)B. Therefore, the following exact intertwining 
relation holds 

V''*^'\t,0)Pf^'\0) = pf^'\t)V''*^'\t,0), k = 0,---,n. (1.8) 

Introducing uj{t) = maxfc=o, -,n ^^k{t) to control the norm of the "dynamical phases" , we prove 
the existence of k > such that for any < t < 1 

U{t,0) = y«*(^)(t,0) + 0(te-'^/Vo"(")'^"/"), (1.9) 

where \\V''*'^^\t,0)\\ = 0(e/(>(«) '^Veg^/^"), for some L> > and < /? < 1. Note that the 
leading term is always meaningful with respect to the exponentially smaller error term. 



3 



In case i3 is a Hilbert space and H{t) is self-adjoint, both the evolution and its approxi- 
mation are unitary and D can be chosen equal to zero. The intertwining identity (|1.8|) and 
show that the transitions between the different subspaces Pj ^^\o)B are exponentially 
small in e, instead of being of order e between the spectral subspaces of H. Constructions 
leading to approximations of this type with exponentially small error term go under 

the name superadiabtic renormalization, according to the terminology coined by Berry in 
this quantum adiabatic context. The first general rigorous construction of this type appears 
in but we shall use that of ||22J. The statement H1.9() is thus very similar to the Adiabatic 
Theorem of quantum mechanics [23 , [101, [Zji P--- and, more precisely, to the subsequent 
exponentially accurate versions in an analytic context provided in j21j . |31j . |22j . j24j . jl6j . 
|17j... or variants therof. However, while the improvement of the error term in (|1.9|) from 0{e) 
to 0{e~'^^^) by considering Pj in place of Pj in the adiabatic context is just that, improve- 
ment, in case there are non-zero nilpotents in the decomposition H1.2|l . it becomes necessary 
to consider Pj and achieve exponential accuracy to get a result. 

This can be understood as follows. As $e(i, 0) commutes with all the P^ ^'^''(0), A; = 0, • • • , n 
for any t we can write 



^i*('){t,0) = Y,Pf^"\0)'^'*^"Ht,0)Pf^'\0) = ^<(")(t,0). (1.10) 

k=0 k=0 



The operator {t, 0) describing the evolution within the fixed subspaces Pj {0)B satisfies 
for j > 1, 

iemf^'\t,0) = {X,{t)pf^'\0) + Dj{t,e) + O{e))<^>f^'\t,0), 
^f^'\0,0)=pf^'\0), (1.11) 

where Dj{t,e) denotes the nilpotent Dj{t,e) = Wi*'^^'>'^ {t)Dj{t)W'i*^^\t). We can write 

^f'-'\t,0) =e-i^*^^(")'^"^f ^'^(t,0), (1.12) 

where the operator ^'^ is essentially generated by a nilpotent. Such adiabatic evolutions 
generated by perturbations of analytic nilpotents are studied in Section IH We show that 
^'^ typically grow when e — > as 

^f^'\t,0) ~ e"/'^', with 0<l3j < 1, (1.13) 

wheras ^f^''\t,0) remains bounded as e ^ iff Dj{t) = 0. The growth in e^^^^ of adiabatic 
evolutions generated by certain nilpotents is already present the works 12 and (2S|- Hence, 
to compensate the exponential growth in 1/e^J of the ^^■'(t, 0)'s which induces transitions 
between the instantaneous eigenspaces of the same order, see the example below, it is necessary 
to push the estimates to exponential order, see (|1.9|) . by trading the Pj's for the Pj This 
requires analyticity of the data, see Section [21 Analyticity is also essential in Section 0] where 
the properties of nilpotent generators and the adiabatic evolutions they generate are studied. 

Let us finally comment on the paper j35j . It addresses, at a theoretical physics level, the 
evolution of master equations describing open quantum systems in which the components of 
the Lindblad generator are slowly varying functions of time. Mathematically, this corresponds 
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to a particular case of problem with a generator containing nilpotents in its decomposition 
(|1.2j) . The authors argue under certain implicit conditions on the evolution, that it is possible 
to approximate U{t, 0) by some operator V^{t, 0) which satisfies the intertwining relation 1)1.81) 
with the instantaneous projectors Pj{t) in place of the approximate projectors Pj ^'^\t) = 
Pj{t) + 0(e). However, as we prove, such a statement cannot be true in general. It does hold, 
however, under the hypotheses of that is when the nilpotent part of the generator in the 
corresponding subspace Pj{t)B is absent, together with an a priori bound on the evolution 
(see also remark iii) at the end of the Section). Or, when the considered spectral subspace 
Pj{t)B is always decoupled from the others, otherwise the error term becomes too large due 
to the growth 1)1.13(1 . An example of this sort is indeed provided in . 

The paper is organized as follows. We close the introduction by the example alluded to 
above and then provide the precise hypotheses and the mathematical statement corresponding 
to our main result. The rest of the paper is devoted to the proof of it. The main steps 
consists in Section^which studies adiabatic evolutions generated by (perturbations of) analytic 
nilpotents. The iterative scheme providing the adiabatic renormalization of |22j is shortly 
recalled in Section [51 The approximations and its properties are presented in Section [HJ 

Acknowledgement: The author would like to thank C. Ogabi, R. Rebolledo and D. 
Spehner for useful discussions. 

1.1 About the effect of nilpotents 

We consider here an explicitely solvable model defined by simple generator with two real val- 
ued distinct eigenvalues possessing a nilpotent in its spectral decomposition. We show that 
this nilpotent induces exponentially increasing transitions (in , /3 < 1) between the instan- 
taneous eigenspaces, thereby underlying the necessity to use superadiabatic renormalization 
to achieve our result. We also identify the approximated projectors Pj that the evolution 
follows. 

Let H he a constant 3x3 matrix in canonical basis {ei, 62, 63} defined by 
/O a 0\ 

H = io (1.14) 
\0 1/ 

and let L be another constant 3x3 matrix defined by 
/ -A;\ 

L = { -k , (1.15) 




where the non-zero scalars a, k will be chosen later on. We set 

S{t) := e"**^, H{t) := S{t)HS-^{t), (1.16) 
and consider the adiabatic evolution U {t, 0) defined for any t E [0, 1] by 

ieU'it, 0) = H(t)U{t, 0), C/(0, 0) = I. (1.17) 
The spectrum of H{t) is {0, 1} and its decomposition reads 

H{t) = S{t){OPo + Do + lPi)S-\t) = OPo(i) + Do{t) + IPi(i). (1.18) 
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where Pq = ei(ei| + 62(62!, -Pi = 63(63! and Dq = a 61(62!- Here {(6j!}j=i_2,3 denotes the 
adjoint basis of {ej}j=i^2,3- 

The operator ft{t) := S^^{t)U{t,0) satisfies 

-it{H-eL)/s 



i£n'{t) = {H- eL)n{t), n{t) = e 
The matrix H — eL is now diagonaUzable and its spectrum is 

{l,+\f£ak,-\f£ak} = {1, A+(£), -A+(£)} = {1, A+(£), A_(£)}, 



(1.19) 



(1.20) 



where denotes any branch of the square root function. The corresponding spectral projec- 
tors are denoted by Pi (£) , P+ {e) and P_ (£) and they are given by 

ek 







l—eak 
21,2 



.0 1 



\+{s)sk 



X+(e)-X-{6) X+(e)-X-(s) (A+(£)-A_(£))(A+(£)-l) 



ek 



A+(£) 



X+(e)-X-{e) X+(e)-X^(e) (A+(£)-A_(£))(A+(£)-l) 















(1.21) 



(1.22) 



and P-(e) has the same expression as P+(e) with indices + and — exchanged. Note that 
P±(e) ~ zba/Veofc as £ ^ 0. Whereas the projectors 



Pl(£) = Pl+0(£) 

Po(£) = P+(£)+P_(£) =Po + 0(£) 

admit expansions in powers of £. Hence 

U{t,0) = 5(t)(e-^*/^Pi(£) +e-^*^+(^)/"P+(£) +e-^*^-(^)/=P_(£)), 
so that, as £ — 0, 



\\u{tM 



-itV eak/e itVeak/e 



2Veak 

which diverges, whatever the nonzero value of ak is. 



(1.23) 
(1.24) 

(1.25) 
(1.26) 



Wc now choose ak < and A-t(e) = ^iy/e\ak\ G fM, for dcfinitcness. Since Pk{t)U{t, O)Pj(O) 
S{t)PkO,{t)Pj, j,k € {0, 1}, where S{t) is independent of e, it is enough to compute Pk^l{t)Pj 
to get the behaviour in e of the transitions between the corresponding instantaneous subspaces. 
We get for t > and Pi = 63(63!, 

/ 

(esl +0(£), as £^0 (1.27) 



Pon{t)Pi 



t-s/\ak\/V^ 



ek 
' 2 





PiQ.{t)PQ = 



(1.28) 



The first formula thus implies that the evolution C/(t, 0) does not follow the instantaneous 
eigenprojcctor Pi(t), whereas the second formula simply reflects the non-generic fact that Pq 
is invariant under H — eL m our example, see the remarks below. 

The model being explicitly solvable, we can readily identify the approximated projectors 
the evolution follows. Setting for j = 0, 1 



P;(i,£) ■.= S{t)P,{e)S-\t) = Pj{t) + 0{e), 



(1.29) 
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we compute by means of (|1.25|) and (|1.23|) 

U{t,0)P*{0,e) = P*{t,e)U{t,0). (1.30) 

Thus the evolution U{t,0) exactly follows the projectors (|1.29)) whereas the transition from 
Pi(0) to Po{t) are exponentially large in 1/y/e. 

Remarks: 

i) If the product ak £ C \ a similar result holds. We took ak < for simplicity. If the 
product ak is positive, the transition does vanish in the limit e ^ 0. This is due to the fact that 
the spectral projector Pq corresponding to the unperturbed eigenvalue of if is of dimension 
2. For the natural generalization of this example with dimPo = d, d > 2, the following holds. 
Generically, the splitting of the unique eigenvalue zero of the nilpotent PqH by a perturbation 
of order e yields d perturbed eigenvalues Xj{e) ~ ae^^'^e^'^^'^^^ , j = 0, • • • , d — 1, a G C, see |26] . 
Hence, one of them has a non vanishing imaginary part that produces exponentially growing 
contributions as e — > 0. 

ii) As already mentioned, this example is non-generic in the sense that Pq is invariant under 
r2(t), see (|1.28p . The choice of non-generic L ()1.15|) was made to keep the formulas simple. 
However, as should be clear from the analysis, a generic choice for L implies an exponential 
increase as e ^ for both Pi^(t)Po and Po^(t)Pi, when ak e C \ M"*". 

iii) The real unperturbed eigenvalues and 1 can be replaced by any different complex numbers 
Ao and Ai without difficulty. The main consequence is that the exponents in H1.25() have to 
be changed according to A-i-(e) i-^ Aq + A-i-(e) and 1 Ai. One can assume without loss that 
Q'Aj < 0, J = 0, 1. Observe that if Aq is real and 9Ai < 0, conclusions similar to H1.27() can be 
drawn. In case SjAq < and Ai is real, the transition Po^}(t)Pi is of order e, e ^ 0. This is a 
case where the results of apply, since the evolution (|1.25)1 becomes uniformly bounded in e 
due to the exponential decay stemming from QAq < 0. 



2 Main Result 

Let us specify here our hypotheses and state our result. 
Let a > and Sa = {z £ C \ dist(z, [0, 1]) < a}. 
HI: 

Let {H{z)}^^g^ be a family of closed operators densely defined on a common domain D of a 
Banach space B and for any ip £ D, the map z i— > H{z)(p is analytic in Sa- 

As a consequence, the resolvent R{z, A) = {H{z) — X)~^ is locally analytic in z for A E 
p{H{z)), where p{H{z)) denotes the resolvent set of H{z). 

H2: 

For t G [0,1], the spectrum of H{t) is of the form a{H{t)) = a(t) U cTo(t), and there exists 
G > such that 

inf dist(cr(t),cro(t)) > G. 

te[o,i] 

Moreover, a{t) = {Ai(t), A2(i), • • • , A„(t)} where Xj{t), j = 1, ■ ■ ■ ,n, n < oo, are eigenvalues of 
constant multiplicity rrij < oo such that 

inf dist(Aj(t), Afc(t)) > G. 
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Let Tj £ p{H{t)) be a loop encircling Xj{t) only. The finite dimensional spectral projectors 
corresponding to the eigenvalues Xj (t) are given by 




The loop Tj can be chosen locally independent of t. It is a classical perturbative fact, see j26| . 
that H2 also holds for the spectrum of H{z) with z £ Sa, provided a is small enough, and 
that the eigenvalues are analytic functions in Sa- By this we mean that the inf^gjo,!] can be 
replaced by inf^g^^ in H2. Hence, (|2.1() also holds for z £ Sa and z Pkiz) is analytic in Sa, 
for A; = 0, . . . , n. Consequently, the eigenilpotents given by Dj{z) = {H{z) — Xj{z))Pj{z) are 
analytic in Sa as well. 

We now state a technical hypothesis needed to deal with evolution operators generated by 
unbounded generators. In case one works with bounded operators only, this hypothesis is not 
necessary. 

H3: 

Let Hoi^t) = Po{t)H{t)Po{t). There exists a complex valued function t i— > Ao(t) such that 
for all t £ [0, 1], HQ{t) + Xo(t) generates a contraction semigroup and G p{{HQ{t) + Ao(t)). 

In other words, H3 says that the solution T(s) = g-^^o{i)s g-iHo{t)s ^j^^ strong equation 
on V idsT{s) = {Ho{t) + Xo{t))T{s) satisfies ||r(s)|| < 1, for all s > 0. By Hille-Yoshida's 
Theorem, H3 is equivalent to the following spectral condition for any t G [0, 1], 

[0,oo) C p{-iHo{t)-iXo{t)) and \\{iHo{t) + iXo{t) + X)-^\\ < \, ^X > 0. (2.2) 

A 

This hypothesis implies that the equation 

iedtUo{t,s)^p = Ho{t)Uoit,s)ip, Uo{s, s)^p = ip, s<t£[0,l], V(/9 G P, (2.3) 

defines a unique strongly continuous two-parameter evolution operator UQ{t,s). It means that 
Uo{t,s) is uniformly bounded, strongly continuous in the triangle < s < t < 1 and satisfies 
the relation Uo{t, s)Uo{s,r) = Uo{t,r) for any 0<r<s<t<l. Moreover, Uo{t,s) maps V 
into T>, also satisfies 

iedsUoit, s)ip = -Uo{t, s)Ho{s)ip, G V, (2.4) 

and is such that HQ{t)Uo{t, s){Hq{s) + Ao)~^ is bounded and continous in the triangle < s < 
t <1. Moreover, see Thm X.70., the following bound holds 

\\Uo{t,s)\\ < e/>^oH^"/^ Vs < t G [0,1]. (2.5) 

Since H{t) = Ho{t) + P{t)H{t)P{t) where P{t)H{t)P{t) is bounded and analytic in t, 
Hypothesis H3 also implies existence and uniqueness of a bona fide evolution operator [/(t, s) 
associated with the equation 

iedtU{t,s)ip = H{t)U{t,s)ip, U{s,s)ip = ip, s<tG[0,l], VipeV, (2.6) 

see EH, Thm 3.6, 3.7 and 3.11. 
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Theorem 2.1 Assume HI, H2 and H3 and consider U{t,0) defined by \2.b]) . For k = 
0,---,n, letPf^'\t) =Pk{t)+0{e) he defined by ^5^, ITTTl) n.ndVi*^^')(t.O) = 
given by ifi.l]) . \t>.l^) and i5.11]) . Define uj{t) = m.a,'Xk=o,---,n^^k{i)- Then, there exists 

a constants k > such that for any < t < 1 

g- uj{u) du/ejj^^^ Q^p?* (e) ^ ^{u) du/eyq*{e) Q)p?' (e) 

+0(te-'^/" sup ||e-^''"(")'^"/"y«'(^)(s,0)Pf ^"^(0)11), 

0<s<t 

with 

Vi*^'\t,0)Pf^'\0) = pf'''\t)V''*^'Ht,0), k = 0,---,n. 
Moreover, for all k > there exists < < 1, > 0, and dj. > 0, with do = 0, such that 

\\V'^*^^\t,0)Pf^'\0)\\ < Cke'^k/e^k^^f'^Xk(.u)du/e^ 

with dj = 0, if and only if Dj{t) = 0, j G {1, • • • , n}, in hl.^) . 
As a direct 

Corollary 2.1 Under the hypotheses of Theorem there exists k > 0, < /3 < 1, and 
D > such that 

where V'i''^''\t,0) = O{e^o ^^^) dn/e^D/ei'y 
Remarks: 

0) The equivalent results hold if the initial time is replaced by any < s < t, mutatis 
mutandis. See Subsection 16.11 

1) As is obvious from the formulation, the natural operators to control are e' 
and e-/o'^(")'^Vey9'(^)(t,0). 

ii) As particular cases of Theorem 12.11 we recover the results of ^2] j IHll ? HH] j HHj ) Q • 

iii) In case k is sufficiently large, the different components of the leading order term have 
amplitudes whose instantaneous exponential decay or growth rates in 1/e may change with 
time. More precisely, assume that 

QXk{u)du> / Lo{u)du — K, VtG[0, 1], and V/c = 0, • • • , n. (2-7) 

JO 

This can be achieved by perturbating weakly a generator for which all Afc are real valued, for 
example. Then, for any initial condition 

n n 

V = ^Pf^'\0)^^^ip,{e)eV, (2.8) 

A:=0 A:=0 

we get 

n 

where the error term is exponentially smaller than the leading terms. Each term of the sum 
decays or grows as e — > with an instantaneous exponential rate given by 9 \].{u)du/ e. 
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Depending on the functions t Q Jq Xk{u)du/e, the index of the component which is the 
most significant may vary with time. 

iv) In case all Afc(t) are real, k = 0, ■ ■ ■ ,n, and H{t) is diagonalizable, we can take di. = for all 
k = 0, ■ ■ ■ ,n, and uj{t) = 0. The evolution U and its approximation are then uniformly 

bounded in e and differ by an error of order e~'^^^. Theorem 12. II thus generalizes Thm 2.4 in 
|23j in the sense that we allow permanently degenerate eigenvalues Xj{t), whereas they were 
assumed to be simple in ^j. 



3 Preliminary Estimates 

We start by recalling a perturbation formula for evolution operators that we will use several 
times in the sequel. 

Let {^(i)}tG[o,i] be a densely defined family of linear operators on a common domain D of 
a Banach space B, and assume t i— > A{t) is strongly continuous. Let B(t) be linear, bounded 
and strongly continuous in t S [0, 1]. Assume there exist two-parameter evolution operators 
r(t, s) and S{t, s) associated with the equations 

idtT{t,s)(p = A{t)T{t,s)(p, T{s,s)=I, (3.1) 
idtS{t,s)^ = iA{t)+B{t))S{t,s)^, S{s,s)=I, (3.2) 

for all (/? G 2? and s < t G [0, !]• Then, for any ip £ D, and any r < s < t £ [0, 1], 

idsiT{t, s)Sis, r)^) = Tit, s)Bis)Sis, r)^, (3.3) 

so that by integration on s between r and t, 

S{t,r)^ = T{t,r)^-i I dsT{t,s)B{s)S{s,r)^. (3.4) 

J r 

Iterating this formula, we deduce the representation 



S{t,r) = J2i-iT dsi ds2--- 



n>Q 

X T{t,si)B{si)T{si,S2)B{s2)---B{sn)T{sn,r). (3.5) 
Further assuming that T(t, s) satisfies the bound 

||T(t,s)|| < Me/X")'^", (3.6) 
for a constant M and a real valued integrable function u i— > w(u), we get from 1)3. 5() 

\\S{t,s)\\ KMe-^ri^^'^^+Mmumdu^ (3_7) 

As a first application of ()3.7() . we get from ()2.5() a first estimate on C/(t, s) that we will improve 
later on 

||C/(t,s)|| < e-''"* (3.8) 
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4 Nilpotent Generators 



For later purposes, we study here the adiabatic evolution generated by an analytic nilpotent, 
in a finite dimensional space. We assume 

Nl: 

For any z E S^, ISK^z) is an analytic nilpotent valued operator in a linear space B of finite 
dimension such that for a fixed integer d > 0, N{z)'^ = 0. 



The detailed analysis of the properties of analytic nilpotent matrices is performed in Section 
5 of the book |Sj. It is shown in particular that such operators have the following structure. 
For any nilpotent N{z) satisfying Nl in 5^, there exists a finite set of points Zq C Sa', with 
a' < a, and, there exists a family of invertible operators {S{z)}^(zs^,\Zo such that for any 

Z e Sa'\ Zq, 

N{z) = S-\z)NS{z) (4.1) 

with S{z) and S^^{z) meromorphic in Sa' and regular in Sa' \ Zq. The set Zq where N{z) is 
not similar to the constant nilpotent N is called the set of weakly splitting points of N{z). At 
these points, the range and kernel of N{z) change. 

We consider Y{t,s), defined as the solution to 

edtY{t,s) = N{t)Y{t,s), Y{s,s) = I, Vs,tG[0,l], (4.2) 

and estimate the way Y(t, s) depends on e, as e ^ 0. Note that we don't need to impose s <t 
since we deal with bounded generators. 

In case N is constant, with N'-^^^ ^ 0, Y{t,s) = e(*~*)^/^ behaves polynomially in 1/e, 
i.e. like {{t — s)/eY~^, as e — > 0. When N{t) is not constant, one may expect that Y{t,s) 
explodes less fast than e^^^., which is the worst behavious as e ^ for bounded generators. In 
such cases, however, Y{t,s) grows typically faster than polynomially in 1/e, as the following 
example shows. For N{z) given by 

N{t) = (!2 "M, (4.3) 



-t ^ 

we get that the solution Y{t,Q) to (|4.2j) reads 

/ cosh(^) _-i=sinh(^) \ 

= \t.oM^) - JsinH^J cosh(-L)^ ^si^M-^) J ' ^^-^^ 

which behaves as e*/^, when e — > 0. The growth is nevertheless slower than exponential 
in 1/e. We show that the characteristic behaviour of Y generated by an analytic nilpotent 
operator is similar. For later purposes, we actually consider generators given by an order e 
perturbation of a nilpotent. 

Proposition 4.1 Suppose the nilpotent N{t) satisfies Nl and let {^(i)}te[o,i] be a family 
of operators on B. Then, there exist c > and < /3 < 1 such that the solution Y{t, s) of 

edtY{t,s) = iN{t)+eA{t))Y{t,s), Y{s,s)=I, Vs,tG[0,l], (4.5) 
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satisfies uniformly in t, s £ [0,1] 

\\Y{t,s)\\<ce'/''. 

Remarks: 

i) Asymptotic expansions as e — > of solutions to such equations are derived in |32]) IHHji in 
the neighbourhood of points which are not weakly splitting points for N(z). 

ii) In case both and 1 are not weakly splitting points, it is possible to take (3 = {d — l)/d, 
which is the optimal exponent, see the example. As we shall not need such improvements, we 
don't give a proof. 

iii) The adiabatic evolution generated by an analytic nilpotent does not have to grow exponen- 
tially fast in as e — > 0. Consider for example (|4.,3j) and (|4.4j) along the imaginary t-axis. 
However, such evolutions cannot be uniformly bounded in e, as the next Lemma shows, under 
slightly stronger conditions. 

iv) It is actually enough to assume t \\A(t)\\ is uniformly bounded on [0, 1]. 

Lemma 4.1 Assume {-^(i)}t6[o,i] ^-^ ^ family of nilpotents and {^(i)}fG[o,i] ^-^ ci family 
of operators onB. Consider Y{t, s) the solution to i4.5\ ). Then 

sup \\Y{t, < c» N{u) = Vs < u < t. 

e>0 



Proof of Proposition 14. H The proof consists in two steps. First we prove the result for 
generators with more structure and then, making use of the results of Section 5 in |S] on the 
detailed structure of analytic nilpotents, we extend it to the general case. 

Lemma 4.2 Assume N{t) = S'^{t)NS{t) where N satisfies N'^ = and where {5'(t)}(g[o,i] 
is a family of invertible operators. Let {^(t)}tg[o,i] a family of operators and set 
B{t) = S{t)A{t)S^'^{t) + S' {t)S'^{t). Then, there exists c > such that the solution Y{t,s) 
of \4-5\j satisfies 

\\Y{t,s)\\ < \\S-HtmS{s)\\^^,erA^^^^^^(-^^^^^^^^^ V.<tG [0,1]. 

Remarks: 

0) The constant c depends on only. 

1) If s > t, the same estimate holds with \\B{u)\\du in the exponent, 
ii) This Lemma also holds in infinite dimension. 

Proof of Lemma 14. 2t Let Z{t,s) = S(t)Y{t, s)S~^{s). This operator satisfies by construc- 
tion 

edtZ{t,s) = {N + eB{t))Z{t,s), Z{s,s)=I, Vs,te[0, 1]. (4.6) 
Let us compare Z{t,s) with 

Zo(^,s)=e^(*-^)/^ s,te [0,1] (4.7) 
by means of (|I15|) . We get 

Z{t,r) = ^l dsi ds2---l dSn 

X Zo{t,si)B{si)Zoisi,S2)B{s2)---B{sn)Zo{sn,r). (4.8) 



n>0' 
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Consider now 

Zs{s) = e'^^-^^', for 6>0. (4.9) 

This operator is such that there exists a c > 0, which depends on N only, such that 

\\Z5{s)\\ < c/6'^~^ Vs > 0, and 0<6<l. (4.10) 

Indeed, on the one hand, we have for s > sq, with sq large enoug h 11^5(5)11 < Ke-^'s'^-^, where 
K is some constant which depends on N only. Maximizing over s > 0, we get e-^'s'^-^ < 
^i-d {d-iY'-\ ^Yie other hand, for ah < s < sq, we have \\Zs{s)\\ < e^ol'^l', so that if 
< <5 < 1, (ICTl holds with c = max(e''oll^li,K((d- l)/e)'^-^). 
Coming back to ()4.8|) in which we make use of the relation 

Zo{t,s) = Zs{{t-s)/e)e'^'-'y', (4.11) 

and (IHU)) . we get 

ft rsi fs„-i 

\\Zit,r)\\ < e'^'-'-^/'Y, dsi ds2--- dsn 

„\n-'r J r J T 



n>0' 



< 



X \\Zs{{t - si)/£)B{si)Zs{{si - S2)/e)B{s2) ■ ■ ■ B{sn)Zs{{sn - r)/e)\\ 
ce^it-r)/s (c J^* \\B{s)\\ ds/6''-'r 



ni 

n>0 



„„(S(t— r)/e , 



The left hand side is independent of 5, which we can chose as 5 = e^^'^, so that we eventually 
get 

||Z(^,r)||<-^e/'(^+^li«(^)li)^^/^'''-^'^^ (4.13) 

from which the result follows. | 

Let us go on with the proof of the Proposition. If Zq n [0, 1] = 0, Lemma (j4.2j) applies and 
Proposition 14. II holds . If not, there exist a finite set of real points {0 < ti < ^2 < • • • < < 1} 
and a finite set of integers {pj}j=i ... „i such that 

ma^{\\S{t)\\,\\S-\t)l\\S'{t)S-\t)\\)=0{l/{t-tjr^), as t^tj. (4.14) 

Since Y is an evolution operator, we can split the integration range in finitely many intervals, 
so that it is enough to control Y{t, s) for s < t £ [v,w] C M where [v, w] contains one singular 
point only. Call this singular point to ^ind the corresponding integer pQ. 

Assume to start with that v < to < w. Let 5 > be small enough and v<s<to<t<w 
so that we can write 

Y{t, s) = Y{t, to + 5)Y{to + 6,to- 5)Y{to - 6, s). (4.15) 

The first and last terms of the right hand side can be estimates by Lemma l4.2| whereas we 
get for the middle term 

\\Y{to + 5,to- 5)11 < e^^^o-^ muHeAiumu_ ^^^^^^ 
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Altogether this yields 

\\Y{t,s)\\ < c^S-\t)\\\\S{to + 6mS-\to-6mS{s)\\/e'^'-'y' (4.17) 

By (|4.14j) . there exists a constant c (that may change from line to line) which in dependent 
of e such that the pre-exponential factors are bounded by Also, since N{t) is and 

A{t) is C70 on [0,1], 

rto+S ft 

\\N{u)+£A{u)\\du<c5 and / \\B{u)\\du < c/S^'' (4.18) 



to-S 



and similarly for \\B{u)\\. Hence, Y{t,s) satisfies the bound 

\\Y{t,s)\\ < ce'^.po.w-i)/d+i)/(52po^2{d-i)/d)_ ^4^^g^ 



Choosing 6 = 6(e) = e'^tpo+i) in order to balance the contributions in the exponent, we get 
with a suitable constant c 

||y(t,s)|| < ce'^/" '™ ^ /e (''o+Dd . (4.20) 
Picking < /3o < 1, we get for yet another constant c 

\\Y{t,s)\\ <ce'/'''\ (4.21) 

A similar analysis yields the same result in case tQ = u or = w. As there are only finitely 
many weakly splitting points to take care of, taking for (3 < 1 the largest of the f3j, for 
j = 1, • • • , m, we get the result. | 

Remarks: 

i) The proof is valid in arbitrary dimension, assuming only 1)4.14^ at a finite number of points. 

ii) The exponents > in 1)4. 14() need not be integers. 

Proof of Lemma 14. H Let Y{t, s) be a solution to 1)4. 5|) and assume N(u) = for all 
s<u<t. Then ||r(t,s)|| < H ^Wl'^", which shows one implication. We prove the reverse 
implication by contradiction. Assume there exists uq £ [s, t] such that the nilpotent N{uo) ^ 
and \\Y{t, s)\\ < c, uniformly as e ^ 0, for all < s < t < 1. We compare Y{t, s) with 

Zo(t,s) = e^("o)(*-^)/^ (4.22) 

and get the following estimate from (|3.4I) and ()3.5p 

\\Zo{t,s)\\ < ce^/.*II^H-^(«o)+£AWN«/e. (4.23) 

By Taylor's formula, there exists a 6 > such that t — s < 5 implies 

ft 

\\N{u) - N{uo) + eA{u)\\ du < c6{5 + e), (4.24) 



for another constant c. Hence, if t — s < 6, with 6 small enough, 

\\Zo{t,s)\\<ce''"/', (4.25) 
for some c. On the other hand, if t — s = 6 and e <^ 6, we have for some c, 

\\Zo{t,s)\\=c{6/ey-\ (4.26) 

Thus, by letting 6 and e tend to zero in such a way that 6'^ e <^ 6, we get a contradiction 
between 1)4. 25() and ()4.26|) . which finishes the proof of the statement. 
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5 Iterative Scheme 

We present here the iterative construction which leads to the construction of (t, s) de- 

veloped in 221, to which we refer the reader for proofs and more details. The first general 
construction of this kind is to be found in 1311. 



Assume HI and H2 with a > small enough so that H2 holds in Sa- 

By perturbation theory in z G Sa, if zq G Sa and Tj £ p{H{zq)), j = 1, • • • ,n are simple 

loops encircling the eigenvalues Aj(zo)) there exists r > such that for any z £ B{zo,r), where 
B{zo,r) is an open disc of radius r centered at zq, Tj G p{H{z)), 

For z £ B{zo,r), we set 

^^■^"^ = / iHiz)-XrU\^P!liz), Poiz) = P^iz), (5.1) 

K\z) = iY,P',\z)P^,{z). (5.2) 



A;=0 

-0 



The operator K is bounded, analytic and we define the closed operator 

H^{z) = H{z)-eK^{z) on V. (5.3) 

For £ small enough, the gap hypothesis H2 holds for all z £ B{zo,r), and we set for e small 
enough 

= "2^/ (HHz)-X)-U\, P^=I-Y^P}{z) (5.4) 

n 

K\z) = z^P,i'(z)P,i(z). (5.5) 

fc=0 

Note that A: = 0, • • • , n, and K'^ are e-dependent and strongly analytic in B{zQ,r). 

We define inductively, for e small enough, the following hierachy of operators for q > 1 

miz) = H{z)-£K'^-\z) (5.6) 

/• " 
P^{z) = -— {m{z)-X)-U\ P',=l-Y,P]{^) (5.7) 

n 

K'^iz) = ^Y,Pl\^)Pm■ (5.8) 

fc=0 

It is proven among other things in j^, see also that the following holds: 

Proposition 5.1 There exists eo > 0, 6 > and g > such that for all q < q*{£) = [g/s] and 
all z £ B{zQ,r), K'^{z) is analytic in Sa, and 

\\K'^(z)-K'^-\z)\\<bq\(-Y (5.9) 



.eg. 

\\K%z)\\<b. (5.10) 
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Remarks: 

i) As a corollary, for 

q = q*{e) = [g/el (5.11) 
we get the exponential estimate 

\\K'i'^^\z) - K'i*^'^-^{z)\\ < (5.12) 

ii) The values of Eq and g which determines the exponential decay above only depend on 

sup \\{H{z)-\)-^\\<oo, 

see [HI for explicit constants. 

iii) Since Sa is compact, at the expense of decreasing the value of a, we can assume that 
proposition 15.11 holds for any z ^ Sa, with uniform constants 5,eo ^-nd b. 

Before we go on, let us recall a few facts from perturbation theory applied to our setting, 
that will be needed in the sequel. 

Assume q < q*{e) and let A e U^^-JI'j C p{H{zq)) and z G B{zQ^r). We can write for 
e < £0 

(i7«(z)-A)-i = {H{z)-eK'i-^{z)-\y^ (5.13) 
= {H{z) - A)-i + e{H{z) - X)-^K'^-\z){H^z) - X)~^ 
= (I - eiHiz) - X)-'K'^-\z)r\H{z) - X)-\ 
Hence, for any j = 1, • • • , n, 

Pj{z) = P,{z) - ^ ^ {H{z) - X)-'K'i~\z){H^{z) - A)-i dX 

= P,{z)-eR]{z) (5.14) 

is analytic in z and the remainder is of order e, together with all its derivatives. Moreover, 
making use of 

{H{z) - Xy^ = {H{z) - Xo)-\l - (A - Xo){H(z) - X)-') (5.15) 
for Aq in p{H{z)), we can write 

H{z)P]{z) = H{z)Pj{z)+eF]{z) (5.16) 
where Pji^) given by 

H{z){H{z) - Xo)-' I (I - (A - X,){H{z) - X)-')K'^-\z){H^ {z) - X)-'^ (5.17) 
The identity 

H{z){H{z) - Xo)-^ = I + Xo{H{z) - X^)-\ (5.18) 

shows that Fj{z) is uniformly bounded as e — > and analytic. 
As a consequence, we have 

Lemma 5.1 Let Fj be defined by \5.17^ . Then 

H'^{z)Pj{z) = H{z)Pj{z)+e{F]{z)-K'i-Hz)P]{z)) (5.19) 

H'^{z)P',{z) = Hoiz)+siF^iz)-K^~\z)P^iz)), (5.20) 

where FSiz) = - ZU^j(') ■ 
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6 The Approximation 

Let q < q*{e) and consider V^, defined as the solution to 

iedtV^t, s)^ = (HHt) + eK%t))V%t, s)^, (6.1) 
(p£D, F«(s,s)=I, 0<s<t<l. 

As H'^ = H - eir^-i we get that 

n 

H\t)+eK^{t) = Ho{t) + ^Pj{t)H{t)Pjit) + e{Ki{t) - K^-\t)) (6.2) 

i=i 

is a bounded, smooth perturbation of HQ{t). The results of guarantee the existence and 

uniqueness of the solution to (|6.1|1 . Moreover, as is well known [221, |27j . further satisfies 

V\t,s)P^{s) = P-i{t)V'^{t,s), VA; = 0,---,n, < s < t < 1. (6.3) 

In order to show by means of (|.S.7j) that V^, with q = q*{e), is a good approximation of U, 

we need to control the behaviour of the norm of as e — > 0. We split into components 
within the spectral subspaces of P^. Set 

n 

V^{t,s) = V^t,s)P',{s) s.t. V^t,s) = Y,y^{t,s). (6.4) 

fc=i 

Since the projectors {Pj^{s)}k=o,-,n have norms uniformy bounded from above and below in 
s £ [0,1] and e > 0, there exists a positive constant 7 such that 

7-1 max \\V^it,s)\\<\\V%t,s)\\<j max \\V^it,s)\\. (6.5) 



k=0,---,n " ' ' " " ' ' " ' k=0,---,n " ^ 



We have. 



Proposition 6.1 There exist constants > 0, k = 0,1, ■ ■ ■ ,n, dj > and < /5j < 1, 

j = 1, ■ ■ ■ ,n such that for all e < £q, and all q < q*{s), 

ll^o'(^'^)ll = ll^''(*'^)^o(^)ll < C7oe^/>o(")'^"/^ (6.6) 
\\V^^t,s)\\ = \\Vi{t,s)Pj{s)\\ < Cje'^j/e^'e'^^s^^(^)<in/e_ (g_7) 

Moreover, ^6. 7| ) holds with dj = if and only if Dj{t) = in 
Proof of Proposition 16. H 

We first consider VQ{t,s), the part of corresponding to the infinite dimensional subspace 
Pq . Because of ()6.3|) . it satisfies for < s < t < 1 and any ip £ D 

iedtVS(t,s)^ = {{H\t) + eK'i{t))P',{t))V^{t,s)^, VS[s,s) = P^is). (6.8) 

Lemma [5.11 shows that the generator of V^{t,s) is equal to HQ{t) plus a smooth bounded 
perturbation of order e. We can thus compare Vg'(i,s) and U{t, s)Pq{s) by means of ()3.7() . 
The fact that the initial condition is Pq{s) instead of the identity simply multiplies the estimate 
by 11^0 (s) II, so that we get 

||F«(t,s)Po«(s)|| < ||Po'(s)||e^i'>«(")'^"/"C^ < e^/^oWdHAc-o, (6.9) 
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where Cq and Cq = CqSUPss[o,i] ||-Po(^)II uniform in e. 

The control of the remaining components is conveniently done by taking advantage of the 
intertwining relation (|6.,3|) as follows. 

Let W be the bounded operator satisfying the equation 

iW'^'{t) = Ki{t)W'^{t), W^iO) = I. (6.10) 

This operator enjoys a certain number of properties. As K'^ is smooth and bounded, the 
solution is given by a convergent Dyson series, and W'^{t) interwines between P^{0) and P^{t)- 
Morerover, and its inverse map D into D, see [21) . Finally, by regular perturbation theory 
and Proposition EH = + 0{e) so that 

sup < oo. (6.11) 

te[o,il 

0<e<l 

Therefore, the bounded operator defined by 

^i{t, s) = Wity^V^it, s)W''{s), < s < t < 1 (6.12) 
satisfies by construction 

['^'^{t,s),P^{0)] = 0, VA; = 0,---,n VO<s<t<l. (6.13) 
We can thus view 

$J(t,s) = $«(t,s)P/(0), j = l,---,n, (6.14) 

<^l{t,s) = 1>^(t,s)Po'(0) (6.15) 

as operators in the finite dimensional Banach spaces Pj{0)B, for j > 1 and in the infinite 
dimensional Banach space Po{0)B. Moreover, thanks to 1)6. 11() . there exists a constant C such 
that, uniformly inO<s<t<l and e > 0, 

C-'\\V,'{t,s)\\ < \\<!>l{t,s)\\<C\\V,\t,s)\\, k = 0,---,n. (6.16) 

The operator ^"jit, s) satisfies for any ip ^ D 

iedt^]{t,s)^ = W'^{t)-^H^{t)V'i{t,s)W'^{s)P]{m]{t,s) 

= P]{Q)W^{t)~^H'i{t)P]{t)W'i{t))P]{Q)^]{t, s)ip 

^ H]m]{t,s)^, (6.17) 

where the generator Hj{t) is bounded, see Lemma 15.11 In a sense, $^(t,s) describes the 
evolution within the spectral subspaces. Let us further compute with = (Pj)'^ and ()5.14|) 

H'^{t)Pj{t) = P]{t){H{t)P,{t)+e{F]{t)-K''-\t))P]{t) (6.18) 
= P]{t){\,{t)P,{t) + D,{t) + e{F]{t) - K'^~\t))P]{t) 
= \,{t)P]{t) + P]{t)D^{t)P]{t) 

+eP]{t){\,{t)R]{t) + F]{t) - K^-\t))P]{t) 
^ P^(t){\,{t)+D,{t))P]{t)+eJ]{t). 

The last term is bounded, analytic in t and of order e. We will deal with it perturbatively. 
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Equations H6.18() suggests to decompose s), j = 1, ■ ■ ■ , n, as 

= e-'i'>^(^)'^"/"$^«(t,s) (6.19) 
where *j(t,s) : Pj{0)l3 Pj{0)B satisfies 

iedt^]it,s) = PjmW^it)-\D,it) + eJ]it))W^{t)Pj{0)^'^^it,s), (6.20) 
^P,s) = Pj{0), 

where, in the leading part of the generator, 

Dj{t) = W'^{tr^Dj{t)W'i{t) (6.21) 

is analytic and nilpotent with Dj{t)'^^ = 0, with rrij = dimPj(t). However, the restriction of 
bj{t) to P^{0)B, P^{0)Dj{t)P^{0), is not nilpotent. Nevertheless, '^]{t,s) satisfies the same 
type of estimates an evolution generated by a perturbed analytic nilponent does: 

Lemma 6.1 Let ^'jit, s) be defined by i6.1fJ\) . for j = 1, • • • , n. Then, there exist < /3j < 1 
and dj > 0, Cj > such that 

\\^%t,s)\\<Cje^^/''\ (6.22) 

Moreover, the estimate holds with dj = if and only if Dj{t) = in il.^) . 

Proof of Lemma I6.lt Equations (|5.14jl and (|l.,'-{jl allow to get rid of the projectors -Pj(O) in 
(16.20(1 up to an error of order e, 

P]{Q)D,{t)P]{Q) = W'^{tr'P]{t)D,{t)P]{t)W'^{t) = b,(t) + eL]{tl (6.23) 

where 

L]{t) = -W'^{tr'{R]{t)D,{t)P,{t) + P,{t)D,{t)R]{t) 

-eR%t)D,{t)R]{t))w\t) (6.24) 

is analytic and of order e". Since W'^{t)'^^ is analytic and uniformly bounded, the nilpotent 
Dj(t) satisfies Nl uniformly in e > 0, and (|6.20() and (|6.24() show that the generator of ^j(i, s) 
satisfies the hypotheses of Pr op osition 14 . H which yields the estimate. The last statement stems 
from Lemma 14. l1 . | 

It remains to gather (|6.16|) . (|6.19|) and Lemma IHTTl to end the proof of Proposition l6.1l | 

6.1 End of the Proof 

Given Proposition 16.11 we can finish the proof of our main statement as follows. 
Applying ()3.4() to U and V^, we get 

U{t,r) = V''{t,r) + i f V''{t,s){K''{s) - K''-^{s))U{s,r)ds. (6.25) 

Jr 

Let t uj{t) be the continuous function defined by 

uj{t) = max 9Afc(t). (6.26) 

k=0,---,n 
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Applying (|6.25|) to Pjiir) and multiplication by e -^r ^(s)ds/£ g^ygg H6.4j) 

t /■* t 

||e-^'-^(")'^"/"([/(t,r) - y'?(t,r))P^^(r)|| < y He" '^("^ '^"/^^'(i^ _ ^<?-i(s))|| 

X (||e-J'>(")'^"/^(C/(s,r) -y5(s,r))Pf(r)|| + He" J'X")'^"/"T4'^(s,r)||) ds. (6.27) 
Proposition 16. II and the definition of ijo{t) yield for any < r < s < 1 

||e-J'>(")'^"/"F^''(s,r)|| < Cfce'^'^/""', (with do = 0). (6.28) 

Further taking q = (?*(e), (|5.12j) . (|6.5|) show the existence of constants I? > and < /t < 5 
such that 

||g-/>(n)d«/eyg*{e)(i^5)(_^g*(s)(5) _ ^5* (e)-! (5)) || < ebC e'^ ' <^ ' ' < Be-'''\Q.20i) 
Hence, we get using < t — s < 1, 

||e-i'>(")'^"/"([/(t,r) - V'i*^'\t,r))Pl*^'\r)\\ (6.30) 

+5e"'^/^ sup ||e-J'>(")'^"/^(?7(s,r)-y5*(^)(s,r))pf (^)(r)||, 

r<s<t 

from which we deduce that if e is so small that Be'''/^ < 1/2, 
sup \\e-^>(^)'i^/^{U{s,r) - V'i*^'\s,r))pf'^'\r)\\ 

<2Be-^/'{t-r) sup He" J'>(^)'^"/"y/^"^(s,r)||. (6.31) 



r<s<t 



r<s<t 



In particular, our main result follows. For e small enough, for any < r < t < 1, and for all 
= 0, • • • , n, 

e-S>i^)du/ejj^^^g~^pq*{e)^^-^ = e" '^("^ '^"/^if* (t, r) (6.32) 

+0{{t-r)e~'''' sup ||e-^X")'^"/^yf ^"^(s,r)||). 

r<s<i 

We chose to estimate the difference U — V^ applied on the projectors, because the norms of 
the different components V^^ vary with k. Of course, ()6.32p also holds with removed 

and VI' in place of vf^'\ 

Making further use of (|6.28j> in the error term of (|6.32j> . we get (lowering the value of 

< K < 5) 

U{t,r) = V'''^'\t,r) + 0{{t - r)e-''l'e^r-^^^)du/e^^ (g^33) 
where y5*(^)(t,r) = 0(e/>(«) '^"/sgi'/e''), for some < /3 < 1, and L> > 0. ■ 
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